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The Cortisol stress response and the molecular programming of the corticoid axis were characterized for the 
first time during early ontogeny in a Mediterranean marine teleost, the European sea bass (Dicentrarchus 
labrax). Sea bass embryos, pre-larvae and larvae at specific points of development were exposed to acute 
stressors and the temporal patterns of Cortisol whole body concentrations and the expression of genes 
involved in corticosteroid biosynthesis, degradation and signaling were determined. Expression of genes 
(grl, gr2, mr, erf) involved into the corticoid response regulation combined with histological data indicated 
that, although a Cortisol stress response is evident for the first time around first feeding, a pattern becomes 
established in larvae at flexion until the formation of all fins. Moreover, mRNA transcript levels of llp- 
hydroxylase and llp-hsd2 showed a strong correlation with the whole body Cortisol concentrations. 
Concluding, our data reveal the presence of an adaptive mechanism in European sea bass at early ontogeny 
enabling to cope with external stressful stimuli and provide a better insight into the onset and regulation of 
the stress response in this species. 



The teleostean hypothalamic-pituitary-interrenal (HPI) axis is a system comparable with the mammalian 
stress axis (hypothalamus-pituitary-adrenal; HP A), as a result of convergent evolution 1,2 and it is of utmost 
importance in stress regulation as well as for the adaptation and/or acclimation of fish to their dynamic 
environment. Stress response includes the primary response, resulting in the rapid increase of circulating cate- 
cholamines and Cortisol, the secondary leading to changes in several haematological and biochemical parameters, 
and the tertiary response that involves alterations at the whole animal and population level 13,4 . The fish's response 
to stressors may be of either an adaptive nature, allowing for homeostatic recovery, or a maladaptive nature 
having adverse effects on survival, growth, immune response, reproductive capabilities, behavior and general 
fitness 13,5,6 . 

Exposure to stress can have a profound impact on the physiology and health of an organism later in life 7,8 . 
Studies carried out in mammals have shown that glucocorticoids play a key role in the programming of brain 
structures that can alter the responsiveness to stress 911 . In fact, it has been shown that exposure to stressors during 
development results in permanent changes in stress coping phenotypes in mammals 9,12 , birds 13 , amphibians 14 , 
and fish 15 . 

In teleostean fish, Cortisol is the principal corticosteroid and plays an important role in a number of physio- 
logical processes including growth, immunoregulation, maintenance of energy balance, and reproduction 2,1618 . 
During HPI axis activation corticotropin-releasing factor (CRF), produced in the hypothalamic preoptic area 
(POA), stimulates the pituitary gland corticotropes to secrete adrenocorticotropic hormone (ACTH), which 
regulates Cortisol synthesis and secretion. In teleosts, Cortisol plays also a vital role in the maintenance of 
hydromineral balance, as fish cannot synthesize aldosterone, and Cortisol carries out this mineralocorticoid 
function 1,19 . Cortisol enters by passive diffusion into the cells where its action is mediated by the 
Glucocorticoid receptor(s) (GR) and the mineralocorticoid receptor (MR) 20 , a class of ligand- activated transcrip- 
tion factors. During larval development, marine teleosts undergo dramatic changes in morphology, growth and 
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metabolism in order to accomplish their metamorphosis into juven- 
ile fish. Throughout this period, Cortisol regulates osmoregulatory 
function 21,22 and is implicated in the metamorphosis from larvae to 
juveniles 23,24 . 

Studies conducted in European sea bass, Dicentrarchus labrax, and 
other species showed the presence of maternal Cortisol in embryos 
and that de novo Cortisol synthesis starts shortly after hatching but a 
significant elevation in whole body Cortisol in response to a stressor 
becomes obvious days to weeks later, depending on the species 25 30 . 

However, our knowledge on the development of the hypothala- 
mic-pituitary- interrenal (HPI) axis of European sea bass and its 
response to stressors during early ontogeny or the molecular 
mechanisms involved is scarce. To this end, we examined the tem- 
poral patterns of Cortisol and genes related to the corticosteroid 
signaling (grl,gr2, mr, erf), corticosteroid synthesis {llji-hydroxylase) 
and Cortisol metabolism or inactivation (lip~hsd2) at various stages 
during early ontogeny in order to assess the ontogenesis of the cor- 
ticosteroid-signaling pathway. Moreover, we subjected European sea 
bass embryos, pre-larvae and larvae to an acute stressor in order to 
determine any differences in the timing or magnitude of the activa- 
tion of the corticosteroid stress axis and the molecular response at 
each developmental point/stage. 

Results 

Ontogeny of the HPI axis in European sea bass. Histology data 
revealed that the brain of sea bass was evident at hatching (Figure 1 a) 



while a morphological differentiation of the pituitary was possible at 
day 5 post hatching (dph) (Figures lb, Id). At this day the first 
thyroid follicles appear and they increase in number and size as 
fish grow (Figure If). Pituitary gland is clearly visible and diffe- 
rentiated at 30 dph when neurohypophysis and adenohypophysis 
are distinct (Figure le, If). Interrenal and chromaffin tissues are 
both located in head kidney adjacent to the cardinal vein. The 
kidney was present at 3 dph and it was characterized by the 
presence of kidney tubules with a distinct morphology (Figure lc). 
However, morphological differentiation of the interrenal and 
chromaffin tissues was difficult at early stages of development 
using routine histology, and only at 28 dph, i.e. around flexion the 
respective endocrine cells were simultaneously evident. 

Temporal patterns of Cortisol content and gene expression at early 
ontogeny. European sea bass embryos had low basal Cortisol content 
(2 ± 0.7 ng g _1 ) that had declined at hatching (0.6 ± 0.3 ng g _1 ) and 
subsequently slightly increased at mouth opening (1.5 ± 0.2 ngg -1 ), 
but the differences were not statistically significant. The first peak (P 
< 0.001) was observed at first feeding (6.8 ± 1.3 ng g~'), after which 
whole-body Cortisol mean concentrations dropped gradually from 
flexion (4.1 ± 1.2 ng g -1 ; P < 0.05) onwards to the formation of all 
fins (1.9 ± 0.5 ng g -1 ) (Figure 2a). All genes assessed in the current 
study were expressed in all developmental stages examined. Tran- 
scripts of grl (Figure 2b) showed a gradual increase throughout early 
development, with lowest mRNA abundance recorded in embryos 




Figure 1 | Histological analysis (a) Sea bass larvae at day 1 post hatching (dph). The brain covers the majority of the head area. Bar: 400 |im. (b) A sea bass 
larva at 5 dph with the pituitary morphologically differentiated (red arrow). Bar: 400 urn, (c) Kidney of a sea bass at 3 dph showing the distinct 
morphology of the kidney tubules (arrows). Bar: 100 urn, (d) higher magnification of picture (b) with the hypothalamus and the pituitary (arrow). Bar: 
100 um. (e) Brain of a 30 dph sea bass with fully differentiated pituitary. Arrow: Adenohypophysis, Arrowhead: Neurohypophysis. Bar: 100 um. (f) Sea 
bass head at 30 dph. Red arrow points to the developed pituitary, black arrows point to thyroid follicles. Bar: 100 um. 
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Figure 2 | Temporal patterns of Cortisol content and gene expression at early ontogeny. Changes in resting (0 h) whole body Cortisol levels and 
mRNA transcript levels of grl, gr2, mr, erf, lip~hsd2 and 1 If! -hydroxylase during the different developmental points/stages (embryos-EM, hatch-HAT, 
mouth opening-MO, first feeding-FF, flexion-FLX, formation of all fins-FINS). Values are means ± standard error (n=3). Means with different letters 
differ significantly from one another (P < 0.05). 
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and highest at the formation of all fins. Expression of gr2 (Figure 2c) 
was higher in embryos than grl and showed a statistically significant 
increase (P < 0.05) at hatching and mouth opening, followed by a 
second increase at first feeding until the formation of all fins (P < 
0.001). The mRNA abundance of mr (Figure 2d) showed a similar 
pattern to grl with minimum levels in embryos, a slight increase at 
hatching and mouth opening and a statistically significant increase 
(P < 0.001) at first feeding until the formation of all fins. Expression 
of erf (Figure 2e) was detected in low levels in embryos, and then a 
bimodal pattern was observed with a statistically significant higher 
level of transcripts at hatching (P < 0.001) and the formation of all 
fins than at mouth opening (P < 0.05). 1 1 [^-hydroxylase mRNA 
(Figure 2f) showed low levels in embryos, hatching and first 
feeding, a statistically significant increase (P < 0.001) in mouth 
opening and the formation of all fins and a peak at flexion. Finally, 
transcript levels of 1 1 fS-hsd2 (Figure 2g) showed a bimodal pattern of 
changes with high level of transcript in embryos and at first feeding 
and lowest at the other examined points or stages of development. 

Ontogeny of the Cortisol stress response and molecular onset of 
genes related to the HPI axis. Figures 3 shows the Cortisol response 
and Figures 4-6 the expression profile of the different stress-related 
genes prior to (0 h) and after (0.5 h, 1 h, 2 h and 24 h) the 
application of the stressor during early ontogeny (embryos, hatch, 
mouth opening, first feeding, flexion and formation of all fins). There 
was no statistically significant effect of stress on Cortisol levels of 
embryos at 0.5 h, 1 h and 2 h post-stress, while there was a 
significant (P < 0.001) increase at 24 h after the application of the 
stressor, i.e. approximately 48 h following fertilization. An identical 
pattern of changes was observed after stress at hatching and mouth 
opening. A statistically significant (P < 0.05) effect of the stressor on 
whole-body Cortisol concentrations soon after the application of the 
acute stressors was observed for the first time at first feeding, where 
whole body Cortisol increased from the 0 h basal values (6.8 ± 1.3 ng 
g" 1 ) to a maximum at 0.5 h(10.8 ± 1.0 ngg~')andl h (9.6 ± 0.5 ng 
g -1 ) to return to basal values at 24 h (3.5 ± 0.5 ng g~')- At flexion, a 
prolonged peak was observed with basal Cortisol values at 0 h (4.1 ± 
1.2 ng g" 1 ) and peak values (P < 0.05) at 0.5 h (17.1 ± 1.7 ng g" 1 ), 
1 h (19.9 ± 6.5 ng g" 1 ) and 2 h (23.4 ± 3.1 ng g" 1 ) post stress. In 
addition, the magnitude of the stress response was higher (P < 0.00 1 ) 
than the respective at first feeding. Finally, at the formation of all fins, 
Cortisol content at 0 h (1.9 ± 0.5 ng g _1 ) sharply increased at 0.5 h 
(17.1 ± 1.6 ng g _1 ) to reach a maximum at 2 h (33.7 ± 2.7 ng g _1 ) 
after stress and return to basal values at 24 h (3.3 ± 1.4 ng g _1 ). The 
magnitude of the stress response was statistically significant higher 
(P < 0.001) than the respective in first feeding and flexion. 
Transcripts of grl (Figure 4a) showed no statistically significant 
changes following exposure to the stressors at hatching and mouth 
opening. However, in first feeding and flexion there was a 1.9- and 
1.8-fold upregulation (P < 0.001) at 24 h, respectively, while at the 
formation of all fins a significant 1.4-fold downregulation (P < 0.05) 
compared to controls was observed at 2 h post-stress. A similar 
pattern of expression was observed in gr2 transcripts (Figure 4b), 
with stable mRNA abundance during hatching and mouth opening 
and an upregulation (P < 0.05) at 24 h at first feeding and flexion. In 
the formation of all fins, gr2 reaches a maximum at 0.5 h after stress 
and then mRNA levels fall down to a minimum at 2 h. Transcript 
levels of mr (Figure 5a) showed no statistically significant changes 
between the different time points, apart at first feeding were a 
decrease is observed from 0.5 h till 2 h after stress. Exposure to 
stressors did not affect erf abundance around hatching, mouth 
opening and first feeding, however, in flexion there was an 
upregulation (P < 0.001) at 2 h and 24 h post-stress. In all fins, erf 
transcripts peaked, as whole-body Cortisol concentrations, at 0.5 h 
post stress followed by a gradual decrease to the basal levels of control 
at 24 h (Figure 5b). 



1 1 ^-hydroxylase expression pattern showed statistically signifi- 
cant temporal differences at all points/stages of development. In 
particular, at hatching there was a significant increase (P < 0.001) 
at 24 h (i.e. hatching 100% completed) post-stress, while at mouth 
opening and first feeding there was a peak at 1 h and 2 h post-stress 
respectively. In flexion, llfl-hydroxylase abundance remained at 
high levels from 0 h till 1 h post stress after which there was a stat- 
istically significant decrease (P < 0.05) at 2 h. In all fins, the pattern 
of changes resembled that of Cortisol, with low transcripts at 0 h, 1.6 
to 1.4-fold upregulation at 0.5 and 1 h followed by a sharp down- 
regulation at 2 h (Figure 6a). 1 1 fi-hsd2 levels did not show significant 
changes around hatching and mouth opening. However, in first 
feeding there was a gradual decrease in mRNA transcripts from high 
values at 0 h to a minimum at 2 h post-stress. In flexion, llfi-hsd2 
abundance was constant except for a sharp increase at 24 h. 
However, if we exclude the values at 24 h, there was a statistically 
significant increase at 1 h post stress, which is no longer "masked" by 
the high values at 24 h. As in the case of llfi-hydroxylase, at the 
formation of all fins the pattern of llji-hsd2 changes are similar to 
that of Cortisol, with peak values at 0.5 and 1 h post-stress followed 
by a drop to the basal values at 2 h and 24 h (Figure 6b). 

Discussion 

The temporal changes of whole-body Cortisol levels during the early 
developmental stages of teleosts show that the initial maternal 
deposit of Cortisol is depleted during embryogenesis and reaches a 
minimum around the time of hatch and then the larva begins to 
synthesize Cortisol de novo, a pattern which is observed in a similar 
way across a number of species. These results are in agreement with 
results obtained in this study as well as in studies in Japanese floun- 
der 16 , tilapia 38 , rainbow trout 25 , Asian sea bass 39 , common carp 33 , 
gilthead sea bream 27 , Atlantic salmon 40 , zebrafish 28 , and in previous 
work in E. sea bass 30 . 

Grl is present in the embryos with very low transcript levels but 
soon after hatch its expression follows a continuous elevation during 
development. Gr2 abundance in embryos is higher than grl, and its 
expression pattern is characterized by an initial elevation in hatch 
and a second in first feeding followed by a relatively steady express- 
ion, thereafter. Previous immunohistochemical and in situ hybrid- 
ization data also verify the presence of gr2 mRNA transcripts and of 
the glucocorticoid receptor early in development with increasing 
expression towards the juvenile stage 41 . In addition we show, for 
the first time, the expression profile of mr, erf, ll\i-hydroxy\ase and 
ll\i-hsd2 during the early stages of development. Mr is present in 
embryos in very low copies and its expression increases as develop- 
ment proceeds, following a pattern similar to that obtained for mr 
expression during embryogenesis in zebrafish 28 . 

Crf expression is detected in embryos in very small quantities; 
however, there is a peak in mRNA transcripts at hatch that decline 
in mouth opening and then gradually increase from first feeding to 
the formation of all fins, indicating a maturation of the HPI axis. The 
expression of 1 1 [3 -hydroxylase is upregulated immediately before the 
rise in larval Cortisol levels which occurs in first feeding, pointing to 
the activation of the steroidogenic pathway around that time, llfi- 
hsd2 mRNA transcripts follow a transient pattern as they appear at 
high levels in embryos, drop in hatch, reach a maximum in first 
feeding and then drop again at the following developmental stages. 
The high amounts in embryos may be associated with the maternal 
Cortisol deposit that needs to be metabolized and the second peak 
which appears in first feeding coincides with the first peak of Cortisol 
during early development of sea bass, and may reflect the immediate 
response of the corticoid system to the sudden accumulation of 
Cortisol. Apart from embryos where the mRNA abundance of llfi- 
hsd2 and llji-hydroxylase is high, the expression patterns of these 
genes are quite similar to the respective patterns obtained from 
zebrafish 28 . 
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Figure 3 | Ontogeny of the Cortisol stress response. The Cortisol response prior to (0 h) and after (0.5 h, 1 h, 2 hand 24 h) the application of the stressor 
during early ontogeny. Values are means ± standard error (n=3). Means with different letters differ significantly from one another (P < 0.05). 



The acute stress challenge tests didn't result in a Cortisol response 
in embryos, hatch and mouth opening stages, apart from a maximum 
at 24 h post stress. However, further research is needed to clarify 
whether this maximum is a result of a delayed stress response or 
reflects differences in the developmental point/stage. Histological 



data showed that the first appearance of a distinct hypothalamo- 
hypophysial-interrenal axis is observed at first feeding, where a peak 
in whole body Cortisol levels was observed at 0.5 h post stress, fol- 
lowed by a protracted decrease until at 24 h when it reached resting 
levels. These results, in accordance with the molecular data, imply 
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that as early as at first feeding sea bass individuals are capable of 
responding to external noxious stimuli. In addition the first peak 
observed in whole-body Cortisol concentrations at first feeding 
reflects the essential role of Cortisol to carbohydrate and protein 
metabolism towards transition to exogenous feeding. As develop- 
ment proceeds, the magnitude and duration of the response is higher 
and a pattern seems to be established from flexion until the formation 
all fins where Cortisol values reach a maximum at 2 h. This is further 
supported by the histological data showing that while the kidney 
tubules with a distinct morphology were present at hatching, a clear 
morphological differentiation of the interrenal and chromaffin tis- 
sues was possible only at 28days post hatch (i.e. at flexion) using 
routine histology. Similar results were found in cichlid fishes where 
the head kidney from 12 to 30 days after fertilization is functionally 
mixed, with the nephron and developing hemopoietic and endocrine 
(chromaffin and interrenal tissue) 42 . These results indicate that even 
at the stage of first feeding fish are capable of a stress-induced stimu- 
lation of Cortisol and that the HPI axis becomes gradually established 
until the development of all fins. This is in accordance with previous 
work conducted in European sea bass 30 , rainbow trout 25 , the yellow 
perch 26 and the zebrafish 28 . However, this is the first time that the 
exact pattern of Cortisol response following exposure to acute hus- 
bandry stressors is revealed at early developmental stages. 

With the aim to shed light on the molecular mechanisms related to 
the onset of the Cortisol stress response, we carried out qPCR experi- 
ments in order to measure the mRNA transcript levels of genes 
related to the HPI axis, grl and gr2, mr and erf, and genes related 
to the biosynthesis and degradation of Cortisol, llfi-hsd2 and 11 
hydroxylase. During HPI axis activation, grl, grl and mr are the 
mediators of the transcriptional effects of circulating Cortisol on 
target tissues and erf produced in the hypothalamic preoptic area 
(POA), stimulates the pituitary corticotropes to secrete adrenocorti- 
cotropic hormone (ACTH) 43,44 , which in turn stimulates synthesis 
and secretion of Cortisol into the circulation 45 . Grl expression levels 
after application of an acute stressor were not altered in any devel- 
opmental stage apart from all fins, where a down-regulation was 
observed at 2 h post stress. The same is the case of the mRNA levels 
oigr2, as the only response to the stress was detected at the stage of all 
fins, where an increase was observed at 0.5 h post stress followed by a 
down-regulation at 2 h. The down-regulation of grl and gr2 is in 
accordance with data from other studies carried out in sea bass 
exposed to very high stocking densities 46 , in coho salmon 47 , in 
Atlantic salmon 6 , in common carp 48 , in the hippocampus of rats 
exposed to increasing corticosteroid levels 49 , in mouse pups with high 
corticosteroid levels due to 24-h maternal deprivation 50 . The statist- 
ically significant increase of the expression levels of both grl and gr2 
observed at 24 h post stress in first feeding may be related, as with the 
case of Cortisol, either to a delayed stress response or to the role of 
Cortisol in metabolism and neural development for the passage of fish 
to exogenous feeding. The mRNA abundance of mr was not altered at 
any of the developmental stages examined apart from first feeding 
where the expression levels decreased from 0.5 h to 2 h post stress. 
However, as this is not repeated in the later stages and especially at all 
fins where the HPI axis is expected to be more mature, this alteration 
may not reflect a response to stress but rather a suppression related to 
the needs of the developmental stage. Cr/expression pattern at all fins 
follows an increase at 1 h and a gradual decrease until 24 h post 
stress, which is in accordance with the pattern observed for Cortisol 
at this stage, where the peak of Cortisol levels is at 2 h post stress. In 
flexion, there is a statistically significant increase at 2 h that remains 
at high values still at 24 h, indicating that at this stage the erf system is 
being established, but the prolonged expression and the delayed 
response compared to the pattern of Cortisol at the same stage, reveals 
that it is not mature yet. 

During fish ontogenesis Cortisol is a critical hormone when 
changes occur at the metabolic demands of the larvae 21,24,27 and it 



is also implicated in neural development and in the induction meta- 
morphosis 16,51 . The inability of the stress response system to respond 
to the stress-elevated Cortisol levels via the erf and the grs at these 
stages until only at all fins, might be of critical importance for the 
survival of the larvae and the normal progress of the development. 
There are very limited data available about the role of gr and mr in 
fish development. The major mineralocorticosteroid in mammals 
and non-mammalian vertebrates is aldosterone. However, in fish, 
deoxycorticosterone (DOC) is considered to be a mr ligand instead 
of aldosterone, as the latter is not detected in fish, but also Cortisol is a 
high-affinity ligand for mr. A recent study carried out in zebrafish 
showed that both GR and MR are present during embryogenesis and 
suggested that gr plays a more important role after hatching in zebra- 
fish, whereas mr is suggested to be important at the earlier stages of 
development, and that after hatching a ligand other than Cortisol, 
perhaps DOC, may be responsible for mr signaling 28 . Studies in mice 
showed that mice lacking a functional gr survive until birth but die 
shortly thereafter due to impaired lung development 52 and that there 
is no abnormal embryonic development detected in a mutant zebra- 
fish line that does not develop corticotropic pituitary cells 53 . Other 
studies demonstrated that knocking down maternal gr leads to devel- 
opmental defects in mesoderm formation in zebrafish 54 and that gr 
signaling is essential for zebrafish muscle development 55 . In the pre- 
sent study, mr expression profiles during the stress response in first 
feeding seems to give mr a more important role at this stage than that 
of grs, which at the later stages of development is inverted. This is in 
accordance with the results obtained from the other studies men- 
tioned above. 

We tested the hypothesis that the molecular events related to the 
appearance of the Cortisol synthesis pathway are tightly linked to the 
enzymes which take part in Cortisol biosynthesis and degradation. 
Therefore, we quantified the temporal expression of 1 1 /^-hydroxylase 
that generates Cortisol from 11-deoxycortisol 56 and llji-hsd2, an 
enzyme that converts the biologically active Cortisol to the inactive 
cortisone. After the acute stress application the transcript levels of 
1 1 ^-hydroxylase appear statistically significant altered at hatch and 
mouth opening but these changes may not reflect a stress response, 
but these changes could represent the role of this enzyme in gonad 
differentiation, as 1 1 ^-hydroxylase appears to be a key transducer in 
the mechanism of sex determination in fish 57 . However, in first feed- 
ing appears a strong -relation between llfi-hydroxylase transcripts 
and Cortisol increase. This pattern continues also in the later stages of 
development, where mRNA expression of llfi-hydroxylase is upre- 
gulated along with Cortisol levels. This is in accordance with data 
obtained for rainbow trout, where mRNA abundance of llfS- 
hydroxylase also increases in response to an acute stressor 58 . In first 
feeding, where the first response is observed, llfi-hsd2 mRNA tran- 
scripts are at high amounts in the larvae and gradually decline to a 
minimum at 2 h post stress; if compared with the Cortisol pattern, it 
becomes clear that the mRNA levels of this enzyme are correlated with 
the amount of Cortisol present. At the following developmental stages 
this becomes more obvious as the maximum mRNA levels of this gene 
are observed at 1 h post stress, just prior to the Cortisol peaks (at 2 h), 
and at 2 h the transcript levels drop again to the resting values. 

The Cortisol stress response and the regulation of genes related to 
the corticoid axis in combination with histological analysis were 
studied for the first time during early ontogeny in European sea bass, 
Dicentrarchus labrax. Sea bass embryos, pre-larvae and larvae were 
exposed to acute stressors and the temporal patterns of Cortisol whole 
body concentrations and the expression of genes involved in corti- 
costeroid biosynthesis (llfi-hydroxylase), degradation (llji-hsd2) 
and signaling {grl, gr2, mr and erf) were determined. Histological 
data showed that although the kidney was present at hatching, mor- 
phological differentiation of the interrenal and chromaffin tissues 
was difficult at early stages of development using routine histology, 
and only around flexion the respective endocrine cells were simulta- 
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Figure 7 | Sampling design. Ontogeny of the neuroendocrine stress response in European sea bass (DPH: Days Post Hatch). An acute stress test was 
applied in all developmental points/stages. 



neously evident. Whole body Cortisol concentrations showed a 
decline from embryos to hatching, remained at low levels at mouth 
opening and peaked at first feeding. In addition, around first feeding 
an effect of stress was evident for the first time on post-stress Cortisol 
concentrations. As development proceeds, a pattern with a higher 
magnitude and longer duration of the Cortisol stress response was 
established from flexion until the formation of all fins. Expression 
data of genes related to the regulation of the corticoid response (grl, 
gr2, mr and erf) indicated that, although a Cortisol stress response is 
evident at first feeding, the HPI axis seems to be fully mature only at 
the stage of all fins. Moreover, the mRNA transcript levels of 11 /?- 
hydroxylase and lip~hsd2 showed a strong correlation with the 
whole body Cortisol concentrations. Overall, the data indicate that 
fish do respond to external noxious stimuli as early as at first feeding 
but the Cortisol stress response becomes fully functional and mature 
until only at the stage of all fins. In addition, several changes occurred 
in early development (embryos, hatching and mouth opening) may 
resemble, apart from an adaptive to stress role of Cortisol, its implica- 
tion in other important aspects of development and metabolism. In 
conclusion, our data reveal the presence of an adaptive mechanism in 
European sea bass at early ontogeny enabling to cope with external 
stressful stimuli and provide a better insight into the onset and regu- 
lation of the stress response in this species. 

Methods 

Animals and husbandry conditions. Batches of fertilized European sea bass eggs 
were obtained from a private fish farm (DIAS S.A.) and transferred to the installations 
of the Institute of Aquaculture, Hellenic Center for Marine Research (Heraklion, 
Crete). Larval rearing was performed applying the pseudogreen-water technique 31 , in 
500 L cylindro-conical tanks, with an initial density of 100 eggs L _1 in which both 
hatching and rearing took place. Tanks were coupled to a biological filter and were 
initially filled with filtered sea water from a deep well. Aeration was provided by means 
of a wooden diffuser located in the tank center at a rate of 150-200 ml min~\ Larvae 
were held during the whole experimental period under mean (±SD) water 
temperature of 18 (±1,6)°C, dissolved oxygen levels of 7.2 ± 0.8 mgl~\ salinity of 36 
and pH of 7.9 ± 0.3. Food was delivered only when inflated swim bladder was 
observed in more than 80% of the population. Exogenous feeding was based on 
rotifers {Brachionus sp.) at 5 individuals ml -1 enriched with proteins and PUFA 
(INVE Aquaculture S.A., Belgium) until 10 days post hatching (dph) while 



phytoplankton (Chlorella sp.) was supplied until 10 dph. Enriched Artemia nauplii 
(Instar II, EG, Artemia Systems S.A., Belgium) were delivered since 10 dph until 
50 dph at 0.5 to 1.0 individual ml -1 . From 30 dph, larvae were offered dry feed 
(PROTON 2-3, INVE Aquaculture S.A., Belgium) using automatic feeders. On day 
50 dph, the type of the dry feed changed to PROTON 3-5 (INVE Aquaculture S.A., 
Belgium). During larval rearing and pre- weaning a sample of 10 larvae was taken daily 
to determine the morphological characteristics and record of total length while 2 
times per week weight measurements were also performed with a sample of 10 
individuals. The trial lasted until individuals completed the formation of their fins on 
45 days post hatch (dph). 

Experimental design. Samples were collected at six different embryonic and larvae 
phases (embryos, hatching, mouth opening, first feeding, flexion and formation of 
all fins; Figure 7/Table 1), prior to and after the application of an acute stress test. 
Different stressors were applied based on the tolerance of larvae at the particular 
developmental stages (Figure 7). In particular, embryos were exposed to 
transportation stress for 8 hours at a density of 50 g L -1 and then to netting and 
air exposure for 1 min until distribution to the incubation tanks. Pre-larvae 
(hatching, mouth opening, first feeding) were exposed to chasing with a net for 
20 sec and high aeration (1,000-1,500 ml min -1 vs. 150-200 ml min -1 ) for 90 sec. 
Larvae (flexion and formation of all fins) were exposed to high aeration (as 
above), chasing with a net for 20 sec, confinement (collection in beakers), and air 
exposure for 5 sec before being transferred to baskets within a 500 L tank. 
Samples for molecular and endocrine (Cortisol) analysis were collected with a net 
at 0, 0.5 h, 1 h, 2 h and 24 h post-stress, flash frozen in liquid N 2 and stored at 
— 80 D C. All experiments were performed in accordance with relevant guidelines 
and regulations. The laboratories of the Hellenic Centre for Marine Research are 
certified and obtained the codes for breeding animals for scientific purposes (EL- 
91-BIO-04). Furthermore all procedures involving the handling and treatment of 
fish used during this study were approved by the HCMR Institutional Animal care 
and use committee following the Three Rs (3Rs, Replacement, Reduction, 
Refinement) guiding principles for more ethical use of animals in testing, in 
accordance to Greek (PD 56/2013) and EU (Directive 63/2010) legislation on the 
care and use of experimental animals. 

Histological analysis. Sea bass larvae were killed with an overdose of anesthetic 
(ethylene glycol monophenyl ether, Merck, 807291) and fixed in buffered formalin. 
Samples were dehydrated in a 70-95% ethanol series and embedded in glycol 
methacrylate resin (Technovit 7100, Heraeus Kulzer, Germany). Serial sections were 
obtained at a thickness of 3-5 um on a microtome (Leica RM2245, Germany) using 
disposable blades. After drying, slides were stained with methylene blue/azure II/basic 
fuchsin 32 and examined under a light microscope in order to record the first 
appearance of the tissues comprising the HPI axis and to describe the relevant organs/ 
tissues. 



Table 1 | Morphological characteristics of European sea bass larvae collected at various stages during early ontogeny 
Development Description DPH Total Length (mm) 



Embryos 
Hatching 
Mouth opening 
First feeding 
Flexion 
Fins 



70% of embryos in 50% epiboly stage 

70% of embryos are hatched 

Mouth opens, complete yolk sac absorption 

First day of exogenous feeding 

65% completed the notochord flexion 

All fins have been developed 



-2 
0 
6 
9 
29 
44 



4.71 
5.25 
11.1 
15.48 



: 0.09 
: 0.19 

:0.44 
0.21 
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Whole body Cortisol. Samples were homogenized according to Stouthart et al. 33 . 
Cortisol was measured in duplicate using a RIA in a 96-well plate according to 
Gorissen et al. 34 . All wells except the 'non- specifics' received 100 u.1 Cortisol antibody 
(Cortisol Antibody[xm210] monoclonal and IgG purified (Abeam) and were 
incubated overnight at 4 C C. Subsequently, the plates were washed three times with 
200 ul/well wash buffer and 100 ul blocking buffer was added to each well in order to 
block the non-specific sites. Plates were covered and incubated for one hour at 37°C. 
After the incubation, 10 ul of standard (4 pg-2048 pg cortisol/10 ul) or 10 ul of 
undiluted homogenate was added to designated wells and 10 ul assay buffer was 
added to the non-specifics and B 0 . All wells received 90 ul (333 Bq) of 3H- 
hydrocortisone (PerkinElmer, USA) solution and plates were incubated at room 
temperature for 4 hours. The plates were then washed three times with wash buffer 
and after the final wash step, all wells received 200 ul of 'Optiphase hisafe-3 
scintillation liquid' (PerkinElmer, USA). Beta-emission was quantified by a 3 min 
count per well using a Microbeta Plus (Wallac/PerkinElmer, USA). 

RNA purification and cDNA synthesis. Samples of embryos, pre-larvae and larvae 
were let to thaw on ice, disrupted and homogenized using the TissueRuptor (Qiagen, 
Hilden, Germany) for 20 s in 600 ul RLT plus buffer (RNeasy Plus Mini Kit Qiagen, 
Valencia, USA). Total RNA was isolated with the RNeasy Plus Mini Kit (Qiagen, 
Valencia, USA). RNA yield and purity was determined by measuring the absorbance 
at 260 and 280 nm using the Nanodrop® ND-1000 UV-Vis spectrophotometer 
(Peqlab, Erlangen, Germany), and its integrity was tested by electrophoresis in 1% 
agarose gels. Reverse transcription (RT) was carried out using 1 ug RNA with 
QuantiTect Reverse transcription kit (Qiagen). 

Primer design. Primers for Glucocorticoid Receptor 1 (grl ), Glucocorticoid Receptor 
2 (gr2), Miner alcorticoid Receptor (mr), Corticotropin Releasing Factor (erf), 
eukaryotic Elongation Factor 1 (eEFla), 40S Ribosomal protein S30 (Fau) and 
ribosomal 18S RNA (18S) were obtained by a previous work of our group 30 ' 35 . Primer 
design for steroid ii/?-hydroxylase (llft-hydroxylase) was based on the available 
sequence with accession number AF449173.2 36 . Forward and reverse primers for 11- 
/?- Hydroxys teroid Dehydrogenase type II (llfi-hsd2) were designed based on the 
conserved regions as revealed by multiple sequence alignments of other teleost fish. 

In the case of 1 1 ^-hydroxylase, the forward primer (ll|3_Fwd) was a 21-mer with 
the sequence 5' -GGAGGAGGATTGCTGAGAACG- 3' and the reverse primer 
( 1 1 P_Rev) an 1 9-mer primer with the sequence 5 ' - AG AGGACGAC ACGCTG AGA- 
3'. For llfi-hsd2, the sequence of the forward primer (hsd_Fwd) was 5' - 
CACCCAGCCACAGCAGGT- 3' and the reverse primer (hsd_Rev) had the sequ- 
nece 5'-ACCAAGCCCCACAGACC- 3'. The products of each primer pair were 
further checked with sequencing in order to confirm that they amplify the desired 
genes. 

Quantitative real-time PCR (qPCR). The mRNA expression of genes encoding for 
grl, gr2, mr, erf, 11 j3 -hydroxylase and ll(3-hsd2 was determined with quantitative 
polymerase chain reaction (qPCR) assays using the KAPA SYBR® FAST qPCR Kit 
(Kapa Biosystems). Reactions were cycled and the resulting fluorescence was detected 
with MJ Mini Thermal Cycler (Bio-Rad) under the following cycling parameters: 
95 C for 3 min (HotStarTaq DNA Polymerase activation step), 94°C for 15 s 
(denaturation step), 60°C for 30 s (annealing step), 72°C for 20 s (extension step), 40 
cycles (step 2-step 4). Levels of grl,gr2, mr, erf, 11 ft -hydroxylase and llfi-hsdl mRNA 
were normalized based on the reference genes 18S, eEFla and Fau. A relative standard 
curve was constructed for each gene, using 4 serial dilutions (1 : 5) of a pool of all 
cDNA samples. We also performed geNORM analysis 37 in order to validate which are 
the most suitable reference genes to serve as an internal control and we concluded to 
eEFla and 18S. 

Statistical analysis. All statistical analyses were performed with SigmaPlot 11.0 
(Jandel Scientific). All data are presented as means ± standard error of the mean 
(SEM). Data were initially screened for normality and homogeneity. Statistical 
comparisons of temporal patterns of Cortisol and gene expression of unstressed 
specimens (0 h) between the different developmental stages were made using one- 
way ANOVA. Statistical comparisons of Cortisol content and gene expression 
between the different time points following exposure to a stressor and the various 
developmental points/stages tested were made using two-way ANOVA. Holm- 
Sidak's honestly significant difference test for multiple comparisons was used to 
determine significant differences among groups. The significant level used was P < 
0.05. 
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